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Abstract — We analyze fading relay networks, where a single- 
antenna source-destination terminal pair communicates through 
a set of half -duplex single-antenna relays using a two-hop proto- 
col with linear processing at the relay level. A family of relaying 
schemes is presented which achieves the entire optimal diversity- 
multiplexing (DM) tradeoff curve. As a byproduct of our analysis, 
it follows that delay diversity and phase-rolling at the relay level 
are optimal with respect to the entire DM-tradeoff curve, provided 
the delays and the modulation frequencies, respectively, are chosen 
appropriately. 

I. Introduction 

Efficiently utilizing the available distributed spatial diversity 
in wireless networks is a challenging problem. In this paper, we 
consider fading relay networks, where a single-antenna source- 
destination terminal pair communicates through a set of K 
half-duplex single-antenna relays. We assume that there is no 
direct link between the source and the destination terminals and 
communication takes place using a two-hop protocol over two 
time slots. The source terminal and the relays do not have any 
channel state information (CSI), and the destination terminal 
knows all channels in the network perfectly. 

Previous work: For setups similar to that described above, 
Laneman and Wornell [1] propose space-time coded cooperative 
diversity protocols achieving full spatial diversity gain (i.e., the 
diversity order equals the number of relay terminals). For the 
setup considered in this paper, Jing and Hassibi [2] analyze 
distributed linear dispersion space-time coding schemes and 
show that a diversity order equal to the number of relay terminals 
can be achieved. In [3], assuming the presence of a direct link 
between source and destination, Azarian et al. show that an exten- 
sion (to the multi -relay case) of a protocol previously introduced 
in [4] is diversity-multiplexing (DM) tradeoff optimal. 

Contributions: In this paper, we are interested in a class of 
simple relaying schemes which is based on linear processing at 
the relay level and hence converts the overall channel between 
the source and the destination terminal into a time, frequency 
or time-frequency selective single-input single-output (SISO) 
channel. This is attractive from an implementation point-of-view, 
as it allows to realize distributed spatial diversity through the 
application of standard forward error correction coding over the 
resulting selective-fading SISO channel. The class of relaying 
schemes analyzed in this paper encompasses phase rolling [5], 
[6] and cyclic delay diversity [7] at the relay level. In [7], it 
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is concluded, through simulations, that a if -relay cyclic delay 
diversity system can achieve a diversity gain of K. In [5], it is 
demonstrated that phase-rolling at the relay level can achieve 
second-order diversity. The contributions in this paper can be 
summarized as follows: 

• We introduce a broad family of relay transmit diversity 
schemes based on linear processing at the relay level. 

• While the (numerical) results in [5], [7] are for the case 
of fixed rate (i.e., the rate does not scale with SNR), we 
provide a sufficient condition on the proposed class of relay 
transmit diversity schemes to be optimal with respect to 
(w.r.t.) the entire DM-tradeoff curve as defined in [8]. The 
tools used to prove DM-tradeoff optimality are a method 
for computing the optimal DM-tradeoff curve in selective- 
fading channels, introduced in [9], and a set of techniques 
described in [3]. 

Notation: The superscripts T H and * stand for transpose, con- 
jugate transpose, and conjugation, respectively. Xi represents 
the zth element of the column vector x, and [X] i . stands for 
the element in the ith row and jth column of the matrix X. 
X o Y denotes the Hadamard product of the matrices X and Y. 
rank(X) stands for the rank of X. Tr(X) and ||X|| F denote the 
trace and the Frobenius norm of X, respectively. Ijv is the NxN 
identity matrix. denotes the all zeros matrix of appropriate size. 
We say that the square matrices X and Y are orthogonal to each 
other if (X, Y) = Tr(XY H ) = 0. All logarithms are to the 
base 2. diag(ai, a 2 , . . . , a^v) denotes the N x N diagonal matrix 
with dj on diagonal entry i. The NxN discrete Fourier transform 
(DFT) matrix F is defined as [F] ln = -Le-^C- 1 )*"- 1 ). 
X r~j £A/*(0, a 2 ) stands for a circularly symmetric complex 
Gaussian random variable (RV) with variance a 2 . Let the positive 
RV X be parametrized by p > 0. The exponential order of X 
in p is defined as v = — . f(p) = g(p) denotes exponential 
equality, in p, of the functions /(•) and g(-), i.e., 

lim I0 f /(P) = lim 

p^oo log p p^oo log p 

The symbols >, <, > and < are defined analogously. 

II. System Model 

Preliminaries: We consider a wireless network with K + 2 
single-antenna terminals, where a source terminal S communi- 
cates with a destination terminal V through a set of K half-duplex 
relay terminals IZi (i = 1,2,..., K). For the sake of simplicity, 
we assume that there is no direct link between S and V. The 



channel^] S — > 72-i, with fading coefficient /j, and TZi — > 2?, 
with fading coefficient /i^, (z = 1, 2, . . . , AT), are i.i.d. CAf(0, 1) 
and remain constant over the time-scale of interest. We define the 
column vectors f = [fx fa ■ ■ ■ /k] T and h = [hx h<z ••• hx} 7 '• 

Communication takes place over two time slots. In the first 
time slot, S transmits N symbols consecutively. The relay 
terminals process the received length- A" sequence using a linear 
transformation as described in the signal model below and 
transmit the result during the second time slot to T>, while 
S remains silent. We assume that S and the relay terminals 
do not have CSI, whereas T> knows f,hi (i — 1,2,..., K) 
perfectly. For simplicity, we assume perfect synchronization 
of the entire network and ignore the impact of shadowing and 
pathloss. Throughout the paper, we assume that N > K. 

Signal model: The vectors x, r,, y € C N represent the 
transmitted signal, received signal at TZi, and received signal 
at V, respectively. The vector r j is given by 



1,2,. 



,K 



(1) 



where p denotes the average signal-to-noise ratio (SNR) (for all 
links) and w< is the iV-dimensional noise vector at TZi, with i.i.d. 
CAf(0, 1) entries. The w., are independent across i as well. The 
transmitted signal x obeys the constraint E{x H x} = AT. 

The relay terminal TZi applies a linear transformation accord- 
ing to G^r,;, where the N x N matrix G^ satisfies G^G^ = 
hltf, scales the result and transmits the signal 



ji JV , ai.ui^a mi, i^^un aiiu uuiiaimiD mi, oignui y G^i . This 

ensures that the per-relay transmit power (per dimension) is given 
by p. We emphasize that enforcing a per-relay transmit power 
of p/K, which leads to a total transmit power (across relays) of 
p, does not change the main statements and conclusions in the 
remainder of the paper. The overall input-output relation reads 



y = 



K 



hifiGiX + z 



(2) 



where the effective noise term z (when conditioned on h) is circu- 
larly symmetric complex Gaussian distributed with E{z|h} = 
andE{zz ff |h} = N I N where N' = (l + ^||h|| 2 V 

Since we will be interested in the mutual information (MI) 
between y and x under the assumption that T> knows all the 
channels in the network perfectly, we can divide (j2ji by y/~N^ to 
obtain the effective input-output relation 



y = 



^i+ P (i + \u 2 )U 



K 
2YE 



hifiGiX + z 



(3) 



where z (when conditioned on h) is a circularly symmetric com- 
plex Gaussian noise vector with E{z|h} = and E{zz H |h} = 
I at. In the remainder of the paper, we shall be interested in 



oo case where 



K 



Vi+pa+iNi 2 ) 



the p 

\A+ll h ll 2 iti 
output relation Q as 

y = \/ioH cff x + z 



1+llhll 



With 



^2 h-ifiGi, we can now rewrite the input- 

(4) 



III. Achieving The Optimal 
Diversity-Multiplexing Tradeoff 

Under the assumptions stated in the previous section, it follows 
that the maximum MI of the effective channel in Q is achieved 
by i.i.d. Gaussian codebooks. The corresponding MI is given by 

1 N ^ 

/(y; x|H cff ) = — ]T log(l + pA„(H eff H« )) (5) 

n=Q 

where the factor 1 /2 is due to the half-duplex constraint. 

The DM-tradeoff realized by a family (one at each SNR p) 
of codebooks C r with rate R — r log p, where r £ [0, 1/2], is 
given by the function 



d{r) = — lim 

p — >oo 



logPJjyr) 
logp 



where P e (p, r) is the error probability obtained through max- 
imum likelihood (ML) decoding. We say that C r operates at 
multiplexing gain r. For a given SNR p, the codebook C r (p) € C r 
contains p 2Nr codewords Xj. 

Next, we compute the optimal DM-tradeoff curve, as defined 
in [8], for the effective channel H c ff and provide a sufficient 
condition on the matrices G^ (i = 1,2,..., K) in conjunction 
with a family of codebooks C r (r e [0, 1/2]) to be DM-tradeoff 
optimal. Following the framework in [8], we define the proba- 
bility of outage at multiplexing gain r and SNR p as 



Po{p,r) =P[J(y;x|H eff ) <r log p] 



(6) 



Directly analyzing |6} is challenging as closed-form expressions 
for the eigenvalue distribution of H e g do not seem to be available. 
However, noting that 



/(y;x|H cfr )</ J (y;x|H off ) 



(7) 



where 



7j(y;x|H cff ) = - log 



iiog(i 



N-l 

-y 

n=0 



^ll H off|lF 



•"^(HcffEfJf) 



(8) 



we can resort to a technique developed in [9] to show that 
the DM-tradeoff corresponding to Jj(y; x|H e ff) equals that 
corresponding to 7(y;x|H c ff). The significance of this result 
lies in the fact that the quantity ||H c g lends itself nicely to 
analytical treatment. 

In the following, we will need the N x K code difference 
matrix defined as 



*(Ax) = [GiAx G 2 Ax • • • G K Ax] 



(9) 



where Ax = x x denotes the code difference vector associated 
with the codewords x, x. Our main result can now be summarized 
as follows. 

Theorem 1: For the half-duplex relay channel in Q, the 
optimal DM-tradeoff curve is given by 



B denotes the link between terminals A and B. 



d(r) = K(l — 2r), r £ [0, 1/2]. (10) 



Let {Gi, G 2 , . . . , Gk } be a set of transformation matrices 
and C r a family of codebooks such that for any codebook 
C r (p) £ C r and any two codewords x,x E C r (p) the condi- 
tion rank(€>(Ax)) = K holds. Then, the ML decoding error 
probability satisfies 

P e (p,r) = p- d W. (11) 
Proof: See Appendix [A] ■ 
Discussion: Theorem[T]shows that the DM-tradeoff properties 
of the half-duplex relay channel in Q are equal to the "coop- 
erative upper bound" (apart from the factor 1/2 loss, which is 
due to the half-duplex constraint) corresponding to a system 
with one transmit and K cooperating receive antennas. Noise 
forwarding at the relay level and the lack of cooperation, hence, 
do not impact the DM-tradeoff behavior, provided the matrices 
Gi and the family of codebooks C r are chosen according to 
the conditions in Theorem 1. Azarian et al. [3], assuming the 
presence of a direct link between source and destination, show 
that extending Protocol I in [4] to the multi-relay case by allowing 
only one relay to transmit in a given time slot yields DM-tradeoff 
optimality w.r.t. the entire DM-tradeoff curve. Our results show, 
however, that DM-tradeoff optimality can be obtained even if all 
relays transmit in all time slots as long as the full-rank condition 
in Theorem 1 is satisfied. Another immediate conclusion that 
can be drawn from Theorem[T|is that cyclic delay diversity [7] 
and phase-rolling [5], [6] at the relay level are optimal w.r.t. the 
entire DM-tradeoff curve, provided the delays, the modulation 
frequencies and the codebooks are chosen appropriately. This 
can be seen as follows. We start by noting that the cyclic delay 
diversity scheme [7] can be cast into our framework by setting 
Gi = Pi where denotes the permutation matrix that, 
when applied to a vector x, cyclically shifts the elements in x 
up by i — 1 positions. With 



(Pi,Pj> = 




(12) 



the condition rank(4>(Ax)) = K takes a particularly simple 
form, namely (FAx) fc ^ for all k e {1, 2, . . . , N}. To see this 
note that rank(*(Ax)) = rank(F*(Ax)) and P; = F H A;F 
where 



e mo] e j0i[l] 



J8i[N-l] 



(13) 



with 6; In] 



Next, we have 



A t = diag 

2irn(i-l) 
JV ' 

rank(F*(Ax)) = rank(S [L 1 2 • • • l K ]) (14) 

where E = diag ((FAx)i, (FAx) 2 , . . . , (FAx)at) and [l;] fe = 

e j6>4fc-i]^ _ i 5 2, , JV, 4 = 1,2, . . . K. As a consequence 

of ( 12 1, the columns of the matrix [L 1 2 ■ • • \k ] are orthogonal 
and hence rank(F<I>(Ax)) = K if X has full rank which is the 
case if (FAx) fe ^ for all k e {1,2, . . . ,N}. 

In the case of phase-rolling [5], [6], we have G^ = ^=Ai. 
Again, the condition rank($(Ax)) = K takes a particularly 
simple form, namely (Ax)t ^ for all k € {1, 2, ... , N}. The 
proof of this statement follows by considering 4>(Ax) directly, 
putting rank(4»(Ax)) into the form of the right-hand side of ( 14 1 



and applying the remaining steps in the argument for the cyclic 
delay diversity case. While the (numerical) results in [7], [5] are 
for the r = case, our analysis reveals optimality of cyclic 
delay diversity and phase-rolling for the entire DM-tradeoff 
curve, provided the codebooks satisfy the full-rank condition 
in Theorem 1 . We finally note that cyclic delay diversity and 
phase-rolling are time-frequency duals of each other in the sense 
that the linear transformation matrices for the two schemes obey 
Gi = ^=FP 4 F ff . 

Relation to approximately universal codes [10]: For the half- 
duplex relay channel investigated in this paper, a family of codes 
C r is DM-tradeoff optimal if 



^min(p) > P 



-2r 



(15) 



where p m i n (p) is the smallest eigenvalue of (<fr(Ax)) ff $(Ax) 
over all Ax = x x with x,x E C r (p). This result follows 
immediately from ( p9| in the proof of Theorem 1 . Based on ( [15} , 
we can conclude (using the same arguments as in Sec. IV. A 



in [9]) that any family of codes C r satisfying ( 15 I will also be 
approximately universal in the sense of [10, Th. 3.1]. 

Relation to code design criteria for point-to-point case: We 
conclude our discussion by pointing out that the conditions of 
Theorem 1 guarantee DM-tradeoff optimality in point-to-point 
multiple-input single-output systems as well. 

IV. Conclusions 

We introduced a family of linear relay processing schemes 
achieving the optimal DM-tradeoff curve of half-duplex relay 
channels. Cyclic delay diversity and phase-rolling were shown 
to be (DM-tradeoff optimal) special cases. Our analysis can 
readily be extended to account for the presence of a direct link 
between the source and the destination terminals. Finally, we 
note that the DM-tradeoff framework seems to be too crude 
to quantify potential performance differences between relay 
transmit diversity schemes with different eigenvalue spread of 
the Gramian matrix of the G ;. 

Appendix A 
Proof of Theorem 1 

We start by noting that an upper bound on the DM-tradeoff 
curve can be obtained by applying the broadcast cut-set bound 
[11] to the described network and evaluating the corresponding 
DM-tradeoff for i.i.d. Gaussian codebooks. It is shown in [11], 
[12] that the broadcast cut amounts to a point-to-point link with 
a single transmit and K (cooperating) receive antennas. Taking 
into account the factor 1 /2 loss due to the half-duplex nature 
of the relay terminals, it follows immediately from the results 
in [8] that the DM-tradeoff curve corresponding to the network 
analyzed in this paper is upper-bounded by 

d(r) < K(l-2r), r e [0,1/2]. 

In the following, we shall show that this upper bound is achiev- 
able, despite the lack of cooperation between the relay terminals, 
provided that, for every r £ [0, 1/2], C r (p) £ C r satisfies 



rank(€>(Ax)) = K for all Ax = x — x with x,x G C r (p). 
We start by noting that 



Po(p,r) > F[J] = P[/ J (y;x|H cff ) < rlogp] 



where 



J = {HefflJjCyjxlHeff) < rlogp} 
is defined as "Jensen outage" event. Since 

K K 



P 

l + ||h|| 2 

where h — h o f and the Gramian 
Tt(G x Gf) • 



h H Kh 



K = — 

N 



Tr(G!Gf) 



Tr(G K Gf) 
Tr(G^Gf) 



we have 



2 log U + P r 



h H Kh 



< rlogp 



In what follows, we write \hi\ 2 — p~ Ui and \fi\ 2 — 
p- y i where and Vi are RVs; the choice of this transfor- 
mation will become clear later. Further, we define the events 
A = {ui,u 2l ■ ■ ■ ,u K ,Vi,v 2 , ■ ■ ■ ,v K \ui > 0,Vi > Vi £ 
{1, 2, . . . , K}} and the complementary event A as the event 
where at least one Ui or Vi is negative. Using the law of total 
probability, we can write 

F[J] = F[A] ¥[J\A] +P[A] ¥[J\A] 

and bound F[J] according to 

¥[A] F[J\A] < F[J] < F[A] F[J\A] + P[A] 
F[A] F[J\A] < F[J] < F[A] F[J\A] (16) 
F[J\A] < F[J] < F[J\A] (17) 



where ( fT6[ ) follows from the definition of the Ui and the Uj, their 
independence and by noting that P [A] decays exponentially fast 



in p. The double inequality ( 17 1 results from lim 



logP[-4] _ 
logp 



0. 



We have thus shown that P[JJ= P[i7|-4J. Next, denoting the 
minimum and maximum eigenvalue of K as A m i n and A max , 
respectively, we get the upper bound 



F[J\A] < 
1 



1 



log 1 



An 



log 1 



log 1 



1 + K 



1 + K" " 

K 



< r logp 
2 ) < rlogp 



(18) 



\ ■ K \ 



< r log p 



and the lower bound 

F[J\A] >~ ' 



log (l + pA max ||h|| 2 ) < rlogp 



(19) 



^log ^l+pAmax^l^H^I 2 ^ < rlogp 

\ log (l + A^^Tp 1 ^-"* ) < rlogp 



where the key steps ( 18 1 and ( 19 1 follow from the Rayleigh-Ritz 

K 

theorem [13] and the fact that 1 < 1+ £ p~ Ui < 1 + K for 

Ui > (i = 1,2,..., K) and p > 1. It can be shown that 
the full-rank condition on €>(Ax) implies rank(K) = K and 
therefore A m ; n > 0. We next define the following events 



B 
U 

c 



tv,i, Ui|max(l - Vi - u t ) > j 

\ „maXj(l — Vi 

AminP 





1 


log ( 




2 




1 


log ( 




2 



< rlogp 



1 + K 



where the max is taken over i — 1,2, 
With these definitions, we arrive at 



, K in all three cases. 



nB] 



where d20b follows from 



P[£] < F[J\A] < F[U] (20) 
c n B] < F[J\A] < F[U nB] + F[U n B] 
DB}< P[J\A] < F[U DB} + P[B] 



^ s*l—Vi—Ui 



1 < 



»=1P 



maxi (1 — Vi —Ui) 



< K. 



P 



Now, we can expand P[U n B] as 

F[U nB]= P [0 < max(l - v t - u t ) < 2r + e a 
1 



K 



\hV\hi\ z < 



1 - Fi v^ 



1 l-2r-£i 



A" 



(l-Fi(Vp)) 



K 



(21) 

K 

(22) 



where e a = log(^ )/log p and F^as) = |Ki(|) with K^-) 
denoting the first-order modified Bessel function of the second 
kind. Further, we have 



=(1-Fx(Vp)) 



(23) 



where for ( |22] > and ( |23| > we used the fact that the CDF of the prod- 
uct of two Rayleigh distributed RVs is given by 1 — 2xKi (2x) for 
x > [14]. In the ensuing discussion, all statements involving 
r hold for r S [0,1/2]. 

Combining (|22| and p3), we get 



P[wn6] + P[B] = (l -f^Vp 11 

^ p -K(\-2r) 



■2r-ei 



K 



(24) 



where the exponential equality in ( 24 1 is proved using a Taylor 
series expansion of 1 — ~F\(l/x) around x — and invok- 
ing asymptotic properties of log ( x ) [15, Eq. (4.1.30-31)]. To 
complete the proof, we establish that P[£ H 6] has the same 
exponential behavior (in p) as ¥[U n B] + P [B\ . Using the same 
arguments as in (|2"T|i-( 24 1, it readily follows that 



P[£n6]=P0< max(l 



Ui) < 2r - e 2 



(l-Fi(v^ 



A" 



A' 



(i-Fi(Vp)) 

= p -K{l-2r) 



where e 2 



and hence 



log(AA max ) 
logp 



. We have thus shown that 



p -K{l-2r) < r[J] < p -K{l-2r) 



P[J]=Pj(p,r) = p- K ^- 2r l 

Since Pj(p,r) < Po(p, r) as a result of and since the outage 
probability is a lower bound to the error probability achieved by 
any code [8], we have 



Pj(p,r) <Po(p,r) <P e (p,r). 



(25) 



Following the approach introduced in [9], we now complete the 
proof of the theorem by identifying a family of codes which has 
P e (p, r) = Pj(p, r) and hence results in a DM-tradeoff curve 
which equals the "Jensen" DM-tradeoff curve derived above. We 
start by writing 

P e (p, r) = V[J] P [error | J] + P [error, j] 
< F[J] +P [error, j] . 

Next, we upper-bound P [error, j\ through the union bound 



error 



t,J] 



< 



„2Nr 



(26) 



where we used the fact that the codebook, C r (p), contains p 2Nr 
codewords and P[x — > x, j\ denotes the maximum pairwise 
error probability (over all codeword pairs and all channels in J) 
for ML decoding. With Ax = x — x, we have 




*(Ax)h||| 



(P)l|h|| 



(27) 



4(1 + Ay 

where p m in(p) denotes the minimum eigenvalue of 



(4>(Ax)) H 4>(Ax) and \27\ follows from applying the 
Rayleigh-Ritz theorem. Substituting ||h|| 2 = ^2f =1 p~~ Vi 
into d27li, we have 



^—Vi—Ui 




K 

4(1 + A ^ F 



(P) 



(28) 



where ( 28 1 follows since the event J requires that 
T,f=i p 1 "" 1 "" 1 > p 2r . Finally, inserting Q into Q , 
we get 



P [error, j] 



< 



„2Nr 



exp 



f'r. 



i(p) 2r 



(29) 



4(1 + if)' 

The proof is complete since $(Ax) has full rank for all Ax 
and for all codebooks in C r and hence p m i n (p) > which 
implies that (29 1 decays exponentially in p for all r £ [0, 1/2]. 
Summarizing our results, we obtain 



P e {p, r) < F[J] + P [error, J] 
<¥[J]+p 2Nr exp 
<V[J} = Pj(P,r) 



Mmin(p) 2r 

'4(1 + A) P 



which combined with ( f25] > yields the desired result. 

References 

[ 1 ] J.N. Laneman and G. W. Wornell, "Distributed space-time-coded protocols 

for exploiting cooperative diversity in wireless networks," IEEE Trans. Inf. 

Theory, vol. 49, no. 10, pp. 2415 - 2425, Oct. 2003. 
[2] Y. Jing and B. Hassibi, "Distributed space-time coding in wireless relay 

networks," IEEE Trans. Wireless Comm., vol. 5, no. 12, pp. 3524-3536, 

Dec. 2006. 

[3] K. Azarian, H. El Gamal, and P. Schniter, "On the achievable diversity- 
multiplexing tradeoff in half-duplex cooperative channels," IEEE Trans. 
Inf. Theory, vol. 51, no. 12, pp. 4152-4172, Dec. 2005. 

[4] R. U. Nabar, H. Bolcskei, and F. W. Kneubiihler, "Fading relay channels: 
Performance limits and space-time signal design," IEEE J. Set Areas 
Comm., vol. 22, no. 6, pp. 1099-1109, Aug. 2004. 

[5] I. Hammerstrom, M. Kuhn, and A. Wittneben, "Cooperative diversity by 
relay phase rotations in block fading environments," in Proc. Fifth IEEE 
Workshop on Signal Processing Advances in Wireless Communications 
(SPAWC), July 2004, pp. 293-297. 

[6] P. Kuppinger, "Transformation of distributed spatial into temporal diversity 
by relay phase rotations," M.Sc. Thesis, Imperial College London, Sept. 
2006. 

[7] S. B. Slimane and A. Osseiran, "Relay communication with delay diversity 
for future communication systems," in Proc. IEEE VTC (Fall), Sept. 2006, 
pp. 1-5. 

[8] L. Zheng and D. N. C. Tse, "Diversity and multiplexing: A fundamental 
tradeoff in multiple-antenna channels," IEEE Trans. Inf. Theory, vol. 49, 
no. 5. pp. 1073-1096, May 2003. 

[9] P. Coronel and H. Bolcskei, "Diversity-multiplexing tradeoff in selective- 
fading MIMO channels," in Proc. IEEE ISIT, Nice, France, June 2007, to 
appear. 

[10] S. Tavildar and P. Viswanath, "Approximately universal codes over slow- 
fading channels," IEEE Trans. Inf. Theory, vol. 52, no. 7, pp. 3233-3258, 
July 2007. 

[11] M. Gastpar and M. Vetterli, "On the capacity of large Gaussian relay 
networks," IEEE Trans. Inf. Theory, vol. 51, no. 3, pp. 765-779, March 
2005. 

[12] H. Bolcskei, R. U. Nabar, O. Oyman, and A. J. Paulraj, "Capacity scaling 

laws in MIMO relay networks," IEEE Trans. Wireless Comm., vol. 5, no. 6, 

pp. 1433-1444, Jun. 2006. 
[13] R. A. Horn and C. R. Johnson, Matrix Analysis. New York, NY: 

Cambridge Press, 1985. 
[14] J. Salo, H. El-Sallabi, and P. Vainikainen, "The distribution of the product 

of independent Rayleigh random variables," IEEE Trans. Ant. and Prop., 

vol. 54, no. 2, pp. 639-643, Feb. 2006. 
[15] M. Abramowitz and I. Stegun, Handbook of Mathematical Functions. New 

York: Dover, 1965. 



